The Kondo screening cloud: what can we learn from perturbation theory? 



> 

m 
o 

o 



o 



Victor Barzykin^ and Ian Affleck^'^ 
Department of Physics^ and Canadian Institute for Advanced Research^, 
University of British Columbia, Vancouver, BC, V6T IZl, Canada 

We analyse the role which the distance scale = vf/Tk plays in the single-impurity Kondo 
problem using renormalization group improved perturbation theory. We derive the scaling functions 
for the local spin susceptibility in various limiting cases. In particular, we demonstrate exactly that 
the non-oscillating part of it should be short-range, i.e., vanish for distances r S> I/^f and show 
explicitly that the interior of the screening cloud does not exhibit weak coupling behavior. 
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Although the Kondo effect has been very thoroughly 
studied for over thirty years [Q , comparatively little the- 
oretical work concentrated on the role of spatial corre- 
lations. Part of the problem is that it is inaccessible to 
the Bethe ansatz and difficult to analyse using Wilson's 
renormalization group. No theoretical and experimental 
consensus has emerged on the question of the length scale 
at which screening of the impurity takes place. On one 
hand, in the scaling language the low energy scale 
implies the presence of the exponentially 
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large length scale = vf/Tk] here vp is the Fermi 
velocity, Tk is the Kondo temperature, J is the Kondo 
coupling, and p is the density of states (per spin) . On the 
other hand, in the Knight shift experiments of Boyce and 
Slichter [Q no such scale was observed. It also appears ex- 
perimentally that alloys with impurity concentration 
n ^ (I/S^k)^, i-e., where the intcr-impurity distances are 
much less than display single- impurity behavior. In 
this letter we attempt to clarify these questions using 
renormalization group (RG) improved perturbation the- 
ory. Related theoretical work includes early perturbative 
calculations [5|, and RG approaches ^,0. Recent theo- 
retical work has also addressed these issues. 

In what follows we consider the standard S'imp = 1/2 
Kondo model, 
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The quantity measured in the Knight shift experiments 
is the local spin susceptibility. 
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where S^^^ = S^^^ 



(1/2) / c;r7/;t(r)(T^V(r) is the total 
spin operator of the impurity and conduction electrons. 
The bulk Pauli contribution, xo ~ p/2 has been sub- 
tracted. Recently a scaling conjecture was made ||], sup- 
ported by numerical results, that in the scaling limit, 
rkp 1, T <C Ep, the spin susceptibility has the follow- 
ing form: 
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where X2kp- smd Xun are universal functions of two scaling 
variables ||l^ . This form follows from the rclativistic one- 
dimensional formulation of the Kondo problem ||ll[] . The 
one-dimensional Hamiltonian in terms of the left-moving 
fields is: 



H = VF 



drxpl{r){id/dr)il^L{r) + i;^ AV'I(0)^Vl(0) • S, 
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The local spin susceptibility x(r, T) in this formalism is a sum of uniform and 2kp parts. 



X™(r,T) ^ [vp/T) < [^[(r)^VL(O +^I(-r-)Y^L(-0]^tot >. 



(5) 



and X2kF is given by the same expression with r re- 
placed by — r in the argument of ^p. Here Stot is the 
total spin in the one-dimensional theory: 
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Renormalizability implies that the functions xa (A 
2kp, un) obey equations of the form: 



XA(T.\,D,rTlvp) = 0, (7) 



where D is the ultra-violet cut-off (the bandwidth), A = 
pj is the dimensionless coupling constant, /3(A) is the 13- 
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function and "/a (A) is the anomalous dimension, which is 
a sum of contributions from the local fermion bilinears 
and from S^^^. In this case, both contributions to 
vanish. S^^^ has zero anomalous dimension because it is 
a conserved operator. The fermion bilinear at r 7^ has 
vanishing anomalous dimension because the interactions 
occur only at the origin; only "boundary operators" re- 
ceive anomalous dimensions in such theories Jill . The 
most general solution of this equation is a general func- 
tion of the effective coupling constant at scale T, Ar and 
of the parameter rT/vp- Swapping the dependence on 
At for dependence on T/Tk and using the fact that x^kp 
is real M we obtain Eq. (0). 
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FIG. 1. Cancellation of the uniform part of the local spin 
susceptibility. 

Consider first Xun- From our lower-order perturbative 
analysis we have found that Xun is zero. This fact is in- 
deed quite general, and can be proven to all orders in 
perturbation theory. We will denote as Xtm,imp(^) the 
part of the uniform spin susceptibility for which the im- 



purity piece S'imp of ^tot is responsible, and the other 
part Xun,ci{'r). [These two contributions are defined as 
zero frequency Fourier transforms. Their sum involves 
a conserved operator, S'j^p, and so Fourier transforming 
was equivalent to dividing by r in Eq. (^.] The fact that 
XMn.imp('') vanishes can very easily be seen. In all orders 
of perturbation theory, one integrates over two free elec- 
tron Green's functions, G, connecting the points (r, t) 
and (0,ri) ( see Fig. la). The result of the r-integration 
is zero: 
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Q sin[7rT(T — ti + ir)] sin[7rT(r2 — t — ir)] 
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= 0, 



because after the change of integration variable, t — > 
exp(i2'!TTT), one encounters contour integration with two 
poles on one side. 

The cancellation of Xim,ci(^) is somewhat less trivial, 
since in addition to the graphs with this type of integra- 
tion, there are other graphs. It can be shown that these 
graphs sum up to a graph where the integration of the 
above type is present (see Fig. lb) and therefore also 
vanish. As a result, Xunif,T) exactly vanishes at finite 
r. This function, however, may include a singular contri- 
bution at r = (or, more correctly, at r of the order the 
ultraviolet cut-off in the one-dimensonal theory, which is 
essentially l/kp). 

We will now analyse the 2kp part of the spin suscepti- 
bility X2kF {f, T) . From our analysis of the perturbation 
theory up to third order, we find that the scaling form 
Eq.(^ is indeed obeyed. The scaling function X'ikp can 
be written: 
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where At is given by: 



At = A + A^ \n{D/T) + \^[ln{D/T) - (1/2) ln(D/T) + constant]. 
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Xe is also given by Eq.(^0|), with T replaced by an- 
other effective scale, E{x) — T/[l — exp(— 47ra::)]e. Eq. 
(10) is universal up to a rescaling of the cut-off, D and 
a change in the constant term. We use this freedom to 
redefine D to simplify our expressions. It is important 
to note that the infrared divergences of perturbation the- 
ory are not cut off at low T by going to small r, as was 
first noticed by Gan 0. It is also neccessary to have a 
high T so that At is small. In the third order, these di- 
vergences are associated with the graph shown in Fig. 2. 
Due to the non-conservation of momentum by the Kondo 
interaction, the bubble on the right gives a logarithmic 
T-dependent factor which is independent of r. Thus, at 



low T, the interior of the screening cloud does not exhibit 
weak coupling behavior. x('", T) decreases exponentially, 
c>c e"^'^^, for r ^ vp/T. In the opposite limit, r ^ vp/T, 
X2kp can be expressed as a polynomial in At and In a;, or 
equivalently At and A^, the effective coupling at scale r. 
(Note that A^ = A^ for r < vp/T.) To third order, the 
local spin susceptibility has the following form: 

(A, + cA^) (1-At) , 

X2kp(.XT,Xr) = ——. ^ — (X<C1), (11) 



[r /nvp) 



4T 



where c is a constant. The T-dependent factor is, to the 
order we work, precisely the total impurity susceptibil- 
ity; Xtt{T). This is the total susceptibility less the bulk 
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Pauli term and its value has been determined accurately 
0. At low T it approaches 1/Tk- <C 1 provided that 
r < ix, even if T < Tk- 
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FIG. 2. Singular third-order graph for x{r,T). 

Boyce and Slichter Q have measured the Knight shift 
from Cu nuclei near the doped Fe impurities, at distances 
up to 5-th nearest neighbor. At these very small distances 
of order of a few lattice spacings, they have found em- 
pirically that the Knight shift obeyed a factorized form, 
x(r,T) « f{r)/(T + Tk), with rapidly oscillating func- 
tion /(r) for a wide range of T extending from well above 
to well below the Kondo temperature. This is essentially 
the same form as Eq. (^) at T 3> Tk and low r. We ob- 
serve from Eq. (^) that the factorization breaks down for 
r > vp/T. Our perturbative approach isn't valid unless 
T ^ Tk, so we can't check factorization at low T. This 
question can certainly be addressed for the overscreened 
large-k Kondo problem, where it is possible to obtain 
reliable low-temperature results from the weak-coupling 
perturbative expansion [l^ . The factorized behavior of 
the local spin susceptibility was also obtained in Rcf . psl . 
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So far, we have only discussed the region r » l/kp- 
It is also interesting to consider the integral of x(^) over 
all space. Although the 2k p oscillations tend to cancel 
at large r, there could be a non-zero contribution from 
r < 1/kp. It is useful to consider various static spin 
susceptibilities: 



Xtt[T) = , Xti(T) 
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< ^imp(^)^imp(0) > dr. 
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It is easy to see that the electron spin polarization in 
the presence of an impurity is determined by the spatial 
integral of Eq.(|), or, equivalently, by Xtt{T) - Xti(T). 
Since 5*^^^ is conserved, these susceptibilities obey the 
RG equation, Eq. , with anomalous dimensions deter- 
mined by the anomalous dimension 7iinp(A) of the opera- 
tor ^i^p. For the three different susceptibilities: 7tt = 0, 
7ti = 7imp, and 
results for /3(A) 
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7ii = 27inip. The low-order perturbative 
and 7imp(A) § are: 
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We should note that the RG equations are only simple 
for the choice of spin correlators Eq.(p^). The equations 
become mixed if written in terms of impurity, impurity- 
electron, and electron parts of the spin susceptibility. 
The general solution of these RG equations has the form: 
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Here $j(AT), ^j{^T) are some scaling functions. From 
our third-order perturbative analysis using Wilson's re- 
sult 1^ for XttiT) we have obtained that the functions 
$j(AT) coincide for all three susceptibilities up to and 
including terms of order A^. When A is small, the scale 
factor in Eq. (|lj) can be easily calculated from perturba- 
tive expansion of /3 and 7i,„p: 
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Thus, at least at high temperatures, where our pertur- 
bative calculation of the scaling functions is valid, the 
integrated electronic susceptibility obeys: 



J x(r,T)dr«-^xtt(T). 



(16) 



This result is mostly given by the electron-impurity cor- 
relator, while the electron-electron piece, 

Xee ^ f^dr {Sli{T)Sim) ^ ^Xtt(T), (17) 



is further suppressed by a power of A. In the scaling limit 
of small bare coupling, A ^ 0, the total polarization of 
the conduction electrons vanishes, at least at high tem- 
perature. The oscillating funtion x(^) integrates to 0, 
and Xtt comes entirely from the impurity- impurity part. 
It should be emphasized, however, that the result is non- 
zero at finite bare coupling A. (A typical experimental 
value of A might be 1/ hi[Ep /Tk) ~ .15.) We conjecture 
that the equality of the scaling functions, <I>j(Ar) defined 
in Eq. (^ holds at aU T, so that Eq. (|l|) is true at all 
T and small bare coupling. In particular, the integrated 
electronic susceptibility then vanishes in the scaling limit 
of zero bare coupling at all T. Precisely this result was 
found at T = from the Bethe ansatz [|l^. However, this 
conjecture is not completely consistent with recent work 
of Lesage et al. |l^ which extrapolates to the isotropic 
Kondo Hamiltonian from an anisotropic model. While 
this may well indicate that our conjecture is wrong, it is 
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also possible that there is a problem with the extrapola- 
tion since we find the susceptibilities to be very singular 
in the isotropic limit. 

X2kF{r,T) is a sum of impurity and electron parts, 
X2fci7,imp and X2feF,ci- While the former obeys the RG 
equation Eq.(|^) with — 7iinp, the latter obeys a more 
complicated RG equation due to operator mixing of Sd 
and Simp- X2fcF, imp has the same A-dependent factor as 
in Eq. (nsl) multiplied by a scaling function: 
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where the scaling function 
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is not the same as in Eq. 
weak bare coupling: 



Thus, we obtain for all r at 
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Hence two different scaling functions are present in the 
experimentally measured Knight shift, and their share 
depends upon the gyromagnetic ratios for the impurity 
and conduction electrons. Unlike for the total spin sus- 
ceptibility, X2kF,ei doesn't vanish in the scaling limit of 
zero bare coupling. However, it is small compared to 
X2kp,imp when the bare coupling and x are both small. 

We may summarize the response of the weak-coupling 
Kondo model to a small magnetic field as follows. The 
impurity spin is much more strongly affected by a weak 
Kondo coupling than is the electron gas. This is con- 
nected with the fact that the free impurity susceptibility 
blows up as T — > 0, whereas the free conduction elec- 
tron susceptibility does not. Thus even a weak Kondo 
coupling drastically affects Xtt at low T, causing the di- 
verging Curie susceptibility to level off at 1 /Tr- . On the 
other hand the affect on the electron gas can only become 
appreciable in a long distance, infrared limit. The effect 
of this, however, is not very dramatic, because of the fac- 
tor of 1/r'^, which arises for purely dimensional reasons. 



At short distances, of ©(I/Zcf), the excess polarization 
of the electron gas produced by the Kondo interaction is 
small. This together with the oscillating nature of the 
long-distance polarization gives rise to a small integrated 
excess polarization of the electron gas of 0{\/Tk)- 
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